Bulk-boundary correspondence in Josephson Junctions 
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We discuss Andreev bound states appearing at the interface between two different supercon- 
ductors characterized by different nontrivial topological numbers such as one-dimensional winding 
numbers and Chern numbers. The one-dimensional winding number characterizes d xy and p x wave 
superconductors. The Chern number characterizes chiral superconductors. The number of interfa- 
cial bound states at the zero-energy is equal to the difference between the topological numbers on 
either sides of the Josephson junction. We discuss also relations between properties of the Andreev 
bound states at the zero-energy and features of Josephson current at low temperature. 

PACS numbers: 74.45.+C, 74.50.+r, 74.25.F-, 74.70.-b 



I. INTRODUCTION 

According to the topological classification of matter—, 
a number of unconventional superconductors have been 
categorized in terms of nontrivial topological numbers 2,3 
such as Zi number, Chern number, and one-dimensional 
winding number. The non-centrosymmetric supercon- 
ductor— £ is topologically nontrivial when the amplitude 
of spin-triplet helical-p wave order parameter is larger 
than that of spin-singlet s wave component 6 -. Such 
superconducting phase is characterized by a topologi- 
cal number Zi = 1. The transport properties of non- 
centrosymmetric superconductors are qualitatively differ- 
ent depending on Z^ number—. The spin-triplet chiral-p 
wave superconductivity in S^RuO^^ is characterized 
by Chern numbers n = ±1—. The Chern numbers here 
is referred to as Thouless-Kohmoto-Nightingale-den Nijs 
(TKNN) number in solid state physics 1 - 1 -. The spin- 
singlet chiral-d wave (n = ±2) superconductivity has 
been suggested in Na^CoO^ -y H 2 Oi2~— , heavy fermionic 
compound er 16 ' 17 , graphenei^, and high-T c superconduc- 
tor o 19 : 20 . Unconventional d xy wave symmetry in high-T c 
superconductors and p x wave symmetry in the polar state 
in 3 He are characterized by the one-dimensional winding 
numbe r 21 i 22 which we call Sato number in this paper. 

The unconventional superconductors have subgap An- 
dreev bound states (ABS) at their surface^ 3 - - — , which 
has been known for some time. Such surface state is 
responsible for unusual low energy transport in high-Tc 
superconductors 2 ^—, chiral-p wave superconductor— ~—. 
In particular in spin-triplet superconductors, the surface 
states attract much attention these days because they are 
recognized as Majorana fermion bound states^* — . The 
proximity effect of spin-triplet superconductors is know 
to be anomalous because of the penetration of the Majo- 
rana bound state into a normal meta l 47 : 48 . 

Today the presence of such surface bound state is ex- 
plained in terms of the bulk-boundary correspondence of 
topological superconductivity. According to the bulk- 
boundary correspondence, the number of the surface 
bound state at the zero-energy would be identical to the 
absolute value of topological number defined in the bulk 



superconductor. In fact, this prediction has been con- 
firmed in a number of theoretical studies. In Josephson 
junctions, the validity of the bulk-boundary correspon- 
dence is not clear, when the two superconductors are 
characterized by Sato numbers. 



In this paper, we discuss the number of zero-energy 
ABS at the interface between two superconductors be- 
long to different topological classification by solving the 
Bogoliubov-de Gennes equation analytically. We first 
study the interfacial states between two superconductors 
belonging to different Sato numbers. Since definition 
of the Sato number requires the presence of the time- 
reversal symmetry (TRS) of the junction, the zero-energy 
ABS appears only when the phase difference across the 
junction (ip) is or it. At ip = or n, we confirm that the 
number of the zero-energy ABS's is equal to the difference 
of Sato numbers in the two superconductors consistently 
with the bulk-boundary correspondence of the topolog- 
ical superconductors. We also show that the Josephson 
current at the zero temperature has large values near 
if = or 7T because of the resonant tunneling through 
ABS at the zero-energy. Next we confirmed that the 
number of zero-energy ABS's appearing at the interface 
between two different chiral superconductors is equal to 
the difference in the TKNN numbers in the two supercon- 
ductors. In contrast to d xy and p x cases, the ABS's at the 
zero-energy do not affect the Josephson current between 
two chiral superconductors at low temperature. We also 
discuss the stability of 7r state at the Josephson junctions 
just below superconducting transition temperature T c . 



This paper is organized as follows. In Sec. II, we dis- 
cuss a theoretical model of Josephson junction consisting 
two topological superconductors. In Sec. Ill, we study 
the interfacial ABS between two superconductors char- 
acterized by different Sato numbers. The number of the 
zero-energy ABS and the Josephson effect are studied for 
two chiral superconductors in Sec. IV. We summarize this 
paper in Sec. V. 



2 



II. MODEL 

Let us consider a Josephson junction consisting of two 
superconductors as shown in Fig. [TJ where the electric 
current flows in the x direction and the junction width 
in the y direction is Lj. We apply the periodic bound- 
ary condition in the y direction and consider the limit of 
Lj — » oo. The Bogoliubov-de Gennes (BdG) Hamilto- 




x 



FIG. 1. (color online). A schematic picture of the Josephson 
junction. 



nian in momentum space reads 
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where &j for j — 1 — 3 are the Pauli matrices, <to is the 
unit matrix in spin space, and fj, is the chemical potential. 
In this paper, we consider following pair potentials A(7), 



in Fig.[T] k x = &f cos^{k y — kp sin 7) is the wavenumber 
on the fermi surface in the x (y) direction, and kF is the 
Fermi wave number. The Sato number is defined for each 
moving angle of a quasiparticle 7 and each Nambu space 
in the presence of TRS. For spin-singlet superconductors, 
the BdG Hamiltonian in Eq. ((T|) is block diagonal in two 
Nambu space: Ml and M2. In Ml, spin of electron- 
like (hole-like) quasiparticle is j" (1). On the other hand 
in M2, spin of electron-like (hole- like) quasiparticle is 1 
(|). In this paper, we assume that d vector in the spin- 
triplet symmetry aligns along the third axis in spin space. 
Under this choice, the BdG Hamiltonian in Eq. (JlJ in the 
spin- triplet cases is also decoupled into Ml and M2. In 
Table. I, we summarized the Sato number VF(7) for d X y 
and p x superconductors. 

In the chiral states, n in Eq. ([3]) must be even inte- 
ger numbers for spin-singlet symmetry, whereas it should 

TABLE I. The correspondence between pairing symmetry and 
the Sato number W . The Sato number can be defined only in 
the presence of the time-reversal symmetry for each direction 
of wave vector on the Fermi surface 7 and for each Nambu 
space. The Sato number of s wave case is always zero (i.e, 
W s ~ 0) because s-wave superconductor is topologically triv- 
ial. For spin-singlet d xy symmetry, the Sato number Wd^ v 
depends also on the Nambu space indicated by Ml and Ml. 
The Sato number for spin-triplet p x wave case W Px is always 
unity for all 7 and the two Nambu space. Here the supercon- 
ducting phase is take to be zero. 
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singlet chiral, 
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(3) 



where A is the amplitude of the pair potential, —tt/2 < 
7 < tt/2 is the moving angle of a quasiparticle as shown 



be odd integers for spin-triplet symmetry. The chiral-p, 
d and / wave superconductors are characterized by the 
TKNN number n = ±1, ±2 and ±3, respectively. The 
TKNN number is defined in the absence of TRS. We note 
that the s wave superconductor is topologically trivial. 
Thus both the Sato number and the TKNN one are al- 
ways zero in the s wave superconductor. 



The energy eigen values of Eq. (J) are E = ±E kj ± with E ky± = + |A ± | 2 , A + = A (7), and A_ = A(tt - 7). 
All the pair potentials in Eq. (J3j) satisfy |A+| = |A_|. In such case, the wave functions in the left and the right 
superconductors in Ml are obtained as^ 
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where j = L (R) indicates the left (right) superconductor and ipj is the macroscopic phase of the superconductor. 
The coefficients A, B, C, and D are the amplitudes of outgoing waves from the interface and a and b are those of 
incoming waves. At the junction interface, we introduce the potential barrier described by Vq6(x). The boundary 
conditions for wave function become 



*i(0,y) =*fl(0,y), 



2m 



dx 



Vr(x, V) 



dx 



x->-0' 



V 9 R (Q,y) = 0. 



(8) 



When we calculate the energy of the interfacial ABS, we put a = b = 0. Since we seek the ABS's for \E\ < \Aj\, 
^>l(x,u) (ty R (x,y)) decays at x — > — oo(oo). The decay length is given by the coherence length £o ~ hvp/i^^) with 
vf being the fermi velocity. By using the boundary conditions in Eq. (|8[), we obtain the relation among A, B, C, and 
D as y[A, B, C, D}* = 0, where [•••]* is the transpose of [• ■ • ] and 3^ is a 4 x 4 matrix calculated from the boundary 
conditions. We define that 3(E) is the determinant of y. It is expressed by 



3(E) =(u\ - vls L+ s* L _)(u 2 R - v 2 R s R+ s R -) 

+ \tn\ 2 {^L V R S *R+ S R- + v l s L+s* L _u 2 R - ve lip u L v L s L+ u R v R s R+ - ve~ lip u L v L s* L ^u R v R s R ~} , 



(9) 



where tp — ipz — <^t? (— vr < ip < it) is the phase difference across the junction and t n = cos7/(cos7 + izo) is the 
normal transmission coefficient of junction with zo = Vo/Hvf- The energy of the ABS's is calculated from a condition 
3(E) = 0. 

When we discuss the Josephson current, we first calculate the reflection coefficients of the junction. By eliminating 
C and D using the boundary conditions in Eq. ©, we obtain a relation between (a, b) and (A, B) as 
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where r^ e and r e h are the Andreev reflection coefficients and r e 
Josephson current is calculated based on a formula^ 
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+uu L v L {s L+ - s L -} {u R - V R S* R+ S R -}] 

where uj n — (2n + 1)ttT is the Matsubara frequency. 



Here we explain the definition of v appearing in Eqs. (|9|l- 
(|13[) . When the two superconductors are in the spin- 
singlet symmetry, v become 1 for the two Nambu space 
M\ and Ml. This is also true when the two supercon- 
ductors are in the spin-triplet symmetry. In these cases, 
the energy of ABS obtained from Eq. (j9)) is degenerate 
in the two Nambu space. As a result, the gives rise a 
factor 2 in the Josephson current in Eq. (jlll) . When one 
superconductor is in the spin-singlet symmetry and the 
other is in the spin-triplet one, we take v = 1 in Afl and 
v = -1 in M2. 

The Josephson current can be decomposed into a series 

of 



uu L v L s L -V R s R+ s R - - vu L v L s L+ u 2 R } 



(12) 



(13) 



I 

When two superconductors preserve TRS, /„ is usually 
zero. The normal transmission probability is defined by 



T 



N 



tt/2 



dj cos7|t„| 



ir/2 



(15) 



J = 



DO 

E 

n=l 



J n sin(ny) + I n cos(rup). 



(14) 



The coefficients J n is, roughly speaking, proportional to 
(Tjv) n for s wave junction. The lowest coupling J\ is 
sensitive to the pairing symmetries of two superconduc- 
tors. For instance, J\ = when one superconductor is 
spin-singlet and the other is spin-triplet. 



III. SATO NUMBER 

When the two superconductors are in the s wave sym- 
metry, we obtain the well known results of the energy of 
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where Ao is the amplitude of pair potential at T = and 
N c is the number of propagating channels at the fermi 
level. The amplitude of the critical current at T = is Jo 
in the s wave junctions 52 . There is no zero-energy ABS 
in the s/s junctions because the s wave superconductor 
is always topologically trivial. 

When s wave superconductor is on the left and d xy 
wave one is on the right, the equation for the energy of 
the Andreev bound states is obtained as 



2Ey / A 2 -E 2 - A 2 \t n \ 2 \9 d \ sin(p = 0. 



(20) 




FIG. 2. (color online). The results for s/d xy junctions. The 
energy of the Andreev bound state £ 3 s /d xy i s plotted as a func- 
tion of 7 at zo = 3 in (a). The current-phase relationship is 
shown in (b). 
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The current-phase relationship (CPR) for T < T c be- 



Here we approximately changes 



E 2 to 



y/ A 2 — E 2 with 9 d = sin(27). This approximation is pos- 
sible because the presence or absence of the zero-energy 
ABS is sensitive only to the sign changing of pair poten- 
tial on the fermi surface 2 ^. The energy of the ABS in TVl 
and that in JV2 arc identical to each other. We find that 



^abs _ e > 



/dx 



and 
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with 
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)sgn[sm(— ) 



sinfas/^J = \t n \ 2 9 d smtp. 



< tt/2, 



(21) 
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(23) 



does not become 



Since -w/2 < a s/dwy ^ c s/d ^ 

zero. At \t n \ — > 0, the two superconductors are separated 
from each other and e s , , =0 for all 7. Such zero- 
energy ABS corresponds to the surface bound state of d xy 
wave superconductor. Namely there is one zero-energy 
surface state for each Nambu space for each 7, which is 
a result of the bulk-boundary correspondence of isolated 
superconductors. At \t n \ 7^ 0, a zero-energy state appears 
for each 7 and each Nambu space only when tp = and tt. 
In Fig. Ufa), we show e s s , d as a function of 7 for several 
choices of tp at 20 = 3. The number of the zero-energy 
ABS at tp = is equal to \W s {-y) - W dxy = 1 for each 
Nambu space. The Sato number can be defined only in 
the presence of TRS in the junction under consideration. 
The zero-energy ABS disappear for tp 7^ because the 
Sato number is not be well defined in the absence of TRS. 
In fact, the Josephson current flows through the junction 
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The lowest order coupling is absent (i.e., Ji = 0) because 
s and d xy wave pairing function is orthogonal to each 
other. At T = 0, the Josephson current becomes 
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s/d x 



cos(y>)sgn[sin(<£>)] 



X Jn 
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cos ~i\t n \ 2 \6 d \ 



2Tn J-k/2 ^/l + \t n \ 2 \9 d sm(tp)\ 



(26) 



The resonant transmission through the zero-energy 
ABS's mainly contributes to the Josephson current at 
T = 0. In fact, the Josephson current becomes large at 
tp = and 7r and jumps as shown in Fig. Hfb), where we 
plot the Josephson current as a function of tp. We fix zq 
at 3, which leads to the normal transmission probability 
T N k, 0.07. At T = 0, the amplitude of the Joseph- 
son current is roughly proportional to Tjv as shown in 
Eq. (I26[) . which is a result of the resonant transmission 
through the ABS at the zero-energy. 

The similar conclusion is obtained when we replace 
the spin-singlet d X y wave superconductor by the spin- 
triplet p x wave one. The energy of ABS is obtained as 
j^abs _ £ « an( j £ c _ These energy are given by 



e «/dx H = Asin 
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Vdx, = Ac0S {—2— ) S § n [ Sm ( — 

5m(a s / Px ) = v\tn\ 2 P smtp, 
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with Op = cos 7. We note that there are four disper- 
sion branches in the ABS because of the contributions 
from two Nambu space. The zero-energy ABS appears 
only when <p = and it. The Josephson current in this 
junction is given by Eq. (f24f with 0<j — » 6 P . 
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FIG. 3. (color online). The results for d xy /p x junctions at 
zo = 3. The energy of the Andreev bound state iejj^ ^ i s 
plotted as a function of 7 in (a). The current-phase relation- 
ship of the Josephson current is shown in (b), where results 
for T = 0.1T C is amplified by 10. 

As shown in the Table. HI d xy wave and p x wave sym- 
metries are classified into the same Sato number for 
< 7 < tt/2 in A/1 and -tt/2 < 7 < in Af2. In 
these cases, ABS's may not appear at the zero-energy. 
On the other hand, two zero-energy ABS are expected 
for < 7 < tt/2 in J\f2 and -tt/2 < 7 < in A/1 because 
of \Wd — W Px | = 2. These prediction can be confirmed 
by considering the equation for energy of the ABS, 



E 1 



\tn\< 



E z 



A 2 9 2 - E 2 JA 2 9 2 - E 2 



+vA 2 9 p 6 d cos(tp)] = 0, 



(30) 



where v = 1 for A/1 and v = — 1 for A"2. By substi- 
tuting E — 0, we find \6<i\ + v9d cos <p — 0. There- 
fore, at (p = 0, the zero-energy ABS appear for < 
7 < tt/2 in M2, and for -tt/2 < 7 < in A/1. To 
have expressions of _E ABS , we apply the approximation 

A 2 \9 p 9 d \ 



^A 2 e 2 - E 2 ^JA 2 ( 
we obtain E ABS = 



£ 2 



and ±e 



E 2 . As a result, 
with 



~dxy/Px 



=A|i, 



l v /|^^|cos((^/2)e(^ d ), (31) 
e^ /px =A|t„| v ^^sin(^/2)e(-^ d ), (32) 



where Q(x) is the step function: = 1 for < x and 
= otherwise. In Fig.^a), we show ±e^ < at z = 3 
and <p = 0.47T. When we consider (p — > 0, there are doubly 
degenerate ABS's at the zero-energy for — tt/2 < 7 < in 
A/1 and for < 7 < 7r/2 in A/"2. The doubly degenerate 
dispersionless ABS's at the zero-energy drastically affect 



the Josephson current which is calculated to be 
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°A ^2T Nj _ 7r/2 
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The expression of the Josephson current is given by 
- sin(2ip) 



(33) 
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A V T I 2T, 
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-N 
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C0S 7 |tn| 4 ^^ 
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(34) 



for T < T c . The lowest coupling vanishes because one su- 
perconductor is spin-singlet and the other is spin-triplet 
(i.e., Ji = in Eq. fig]) ). At T = 0, we find 



J d xy / Px =[sin(p/2) - cos( (y 5/2)sgn(v?) 
1 



X Jn 



2T, 



d 1C o S1 \t n \J\8 d 6p\, (35) 



In Fig. EJb), we plot the Josephson current as a function 
of (p at zq — 3. At very low temperature, the Josephson 
current shows unusual current-phase relationship and has 
large amplitude proportional to \/Tn around <p = 0. At 

<£> = 7T, ±e 



dxy/p 



describes the doubly degenerate disper- 
sionless ABS at the zero-energy which also contribute to 
the large Josephson current as shown in Fig. (3Jb). Tun- 
ing of ip at 7T is equivalent to changes the sign of W dxy in 
Table. I. In Fig. [U we show the critical Josephson cur- 
rent as a function of temperature for < T < 0.2T C 
with zq = 3. In this temperature range, the maximum 



value of J s / S is saturates at Jn 



The results for J, 



and Jd^y/px are much smaller than Jo for T > 0.1T C . 
They increases rapidly with decreasing temperature for 
T < 0.1T C . Such effect is known as the low-temperature 
anomaly of the Josephson curren t 28 i 32 . 



IV. TKNN NUMBER 

In this section, we consider a junction which consists of 
two chiral superconductors, where the TKNN number in 
the left superconductor is n and that in the right one is 
to. We can calculate the energy of ABS and the Joseph- 
son current without any further approximations because 
the pair potential is independent of 7 and TKNN num- 
bers. We first show the results of the energy of ABS 
and the Josephson current. Secondly we count the num- 
ber of zero-energy ABS of the junction. Then we discuss 
properties of the Josephson current. 

We show the results for a junction which consists of 
two spin-singlet chiral superconductors. In this case, the 
two TKNN numbers (n in the left superconductors and 
to in the right one) are even integer numbers. The energy 



6 




T/Tr 



of the junction and the spin-triplet chiral superconductor 
occupies the right hand side. The TKNN numbers in the 
left superconductor n is even integers and that in the 
right one m is odd integer numbers. We find the energy 



of ABS is E ABS = 4° with 



: T Asgn[cos(X|°)] sm(X e ± °), 
<°±{(n-m)7-7r/2} 



(44) 
(45) 



cosc£° =(1 - \t n \ 2 ) sin {(n + m)j} + \t n \ 2 vam(p. (46) 
The Josephson current becomes 
A 



T - 1 f j cos(7)|i„| 2 i/cos(^ 

<4o =JOT- > ^7" / 07 -7— To 

A ° v ~ii 2Tn J-*/* 4 sin a™ 



x [sin(2X™)A( e ;°) + sin(2X!°)A(e!°)] . (47) 
where v = ±1 indicates two Nambu space. 



FIG. 4. (color online). The maximum value of the Joseph- 
son current is plotted as a function of temperature. Here we 
choose zq — 3. 



of ABS's are obtained as E ABS = e| e with, 

4 e = ± Asgn[sin(A4 e )] cos(X^), 
= a ee ±(n-m)j 



(36) 
(37) 



cosa ee =(1 - |i„| 2 ) cos {(77 + 777)7} + \t„{ 2 cos(p, (38) 

where < a ee < ir. The Josephson current J ee are also 
expressed in terms of the energy of ABS 
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A 2Tat 



" /2 ^cos(7)|t„| 2 sin(p 



tt/2 



2 sin a ee 



x [sin(2A^ e ),4(4 e ) + sm(2Xt e )A(e e J)] 



(39) 



Next we show the results for a junction where the two 
superconductors belong to the spin-triplet chiral states. 
The two TKNN numbers (n in the left superconductors 
and m in the right one) are odd integer numbers. We 
find the energy of ABS as E ABS = 4° with, 



4° = + Asgn[cos(A£°)] sin(X£°), 
a°° ± (n — 771)7 



A4 = 



(40) 
(41) 



cosa°° =(1 - |rj„| 2 )cos{(77 + 777)7} - |t«| 2 cos^. (42) 
The Josephson current is calculated as 
A 1 



Joo —Jo 



" /2 cos(7)|i n | 2 sin^ 

«7 

., /2 2sma°° 



A 2T N 

x [sm(2Xl°)A(e™) + sm(2X° )A(e°°)] . (43) 

Finally we show the results for a junction where the spin- 
singlet chiral superconductor occupies the left hand side 




n=i m = -3 



n=i m = 3 



FIG. 5. (color online). The left- and the right-hand side of 
Eq. (|50p are plotted with a solid line and two broken lines, 
respectively. We choose n = 1, m = —3, zq = 0.5, and 
tp = 0.2-7T in (a). The parameters in (b) are n = 1, m = 3, 
zo = 3, and tp = 1.37T. 

From the expression of the bound state energy, it is 
easy to count the number of ABS at the zero-energy. For 
example, we consider the junction with both n and m 
being odd integers here. At \t n \ = 0, the energy of the 
energy of the bound state becomes 



e°° = — Asgn(cosn7) sin 717, 
e°° =Asgn(cos 777,7) sin 7777, 



(48) 
(49) 



because of a°° = (77 + 777)7. They represents the disper- 
sion of the surface bound states when the two supercon- 
ductors are separated from each other. It is easy to con- 
firmed that the number of the zero-energy surface states 



7 



of e°° is \n\ and that of e°° is |m|, which is a results of the 
bulk-boundary correspondence for an isolated supercon- 
ductor. At finite |t„|, the solutions of e ao — in Eq. (j40|) 
requires the relation 



tan 



(n — 771)7 



± tan 



(50) 



The left hand side of Eq. ((50)) goes positive infinity 
\n — m\/2 times and goes negative infinity \n — m|/2 
times within the interval of —n/2 < 7 < tt/2. Since 
< a 00 < 7r, tan(% a J on the right hand side remains 
positive value for all 7. Thus the number of the solutions 
of Eq. (|50[) is \n— m\ which corresponds to the difference 
in the TKNN numbers between the two superconductors. 
In Fig. [SJ we plot the left hand side of Eq. (fBT)]) with a 
solid line as a function of 7 and the right hand side with 
two broken lines. In Fig.[5](a), we choose n — 1, m = —3, 
z = 0.5, and ip = 0.2ir. The solid line and two bro- 
ken lines cross four times, which means four zero-energy 
ABS's (i.e., 4 = |l-(-3)|). The number of the solutions 
independent of junction parameters such as z and ip. In 
the lower panel in Fig. [Sfa), we illustrate the schematic 
picture of the chiral edge states at \t n \ = 0. The arrows 
indicate the direction of chiral current. In Fig. [5{b), we 
choose n = 1, m = 3, zq = 3, and tp — 1.37T. The results 
show that there are only two solutions (i.e., 2 = |1 — 3|). 
In the lower panel, there are four chiral edges at \t„ \ = 0. 
The chiral current in the left superconductor flows the 
opposite direction to that in the right superconductor in 
this case. For finite |t„|, two chiral edge current oppo- 
site directions are cancel each other out and two ABS 
remain at the zero-energy. Therefore we confirmed that 
the number of the ABS at the zero-energy is \n — m\ for 
each Nambu space when the TKNN number of two su- 
perconductors are n and m. The similar argument can 
be applied to Eqs. ([36]) and (|44|). The presence of \n — m\ 
zeros has been suggested by the index theorem^. Here 
we count the number of ZES by solving Eq. © explicitly. 

The Josephsson current in Eqs.(|39|) and (|43t have a 
common expression near T c , 



J=Jo- 



A A 



1 



■ sin tp- 



-l 



(n— m) 



(*o), 



(51) 



A T r 4T N 

,7T/2 

In-m(zo) = / d7Cos(7)|i„| 2 cos{(n~ 771)7} (52) 

J -tt/2 

where n and m are both even integer numbers or both 
odd ones. The junction is stable at 7r-state in the case of 
In-m < 0- At zq = 0, the integral for highly transparent 
junctions becomes 



I n -m(ZQ = 0) 



:(-!)" 



(53) 



1 — (n — m) 2 

On the other hand in the tunneling limit zq ^> 1, we find 

(-1)^ 



In-m(z > 1) 



z 2 {(n - m) 2 - 1} {(n - m) 2 - 9} ' 



(54) 



In Table. |TT1 we indicate the stable state of the Joseph- 
son junctions near T c . At n — m = 4, 7r-state is stable 
in the highly transparent limit and 0-state is stable in 
the tunneling limit. Therefore the junction undergoes 
the transition from 7r-state to 0-state when we decrease 
the transparency of the junction. Such junction can be 
realized with chiral-d wave superconductor. 



TABLE II. The stable states of the Josephson junction near T c 
are indicated by '0' or V. Here n and m are both even integer 
numbers or both odd integer numbers. As a consequence, 
n — m becomes an even integer. 



n 


— m 





2 


4 


6 


8 


10 




= 








TV 





7T 





zo 


> 1 











7T 





7T 



In contrast to Sec. Ill, the ABS's at the zero-energy 
do not affect the Josephson current at low temperature. 
At T = 0, Eq. j39) becomes 



Jee —Jo 



SUl ip 

sin(2A 



" /2 cos( 7 )|< n 

-tt/2 



' sintp 



2 sin a ee 
sin(2Af_ e 



|cos(Xi e )| |cos(Ai e |) 



(55) 



Mathematically speaking, the zeros in the |cos(A|. e )| is 
removed. Instead of the zeros of E ABS , minima of sina ee 
in the denominator determine the amplitude of Joseph- 
son current at low temperature. The Josephson cur- 
rent shows the logarithmic dependence of temperature 
at intermediate temperature region between T = and 
T = r e i 35 i 37 . This argument can be applied also to the 
Josephson current in Eq. (|43|) . 

When n is an even integer and m is an odd integer, 
the Josephsson current in Eqs. (1471) near T c becomes 



J eo = - sin(2^)J ^ ( — 



1 



2Ta 



N J -tt/2 



71 72 j cos(7)|i„| 4 cos{2(n-m)7} 
07 



90 



(56) 



where n and m are an even and an odd integer number, 
respectively. The coefficient proportional to — sm(2p) is 
positive in the limit of both zq = and zq 3> 1. 

In Fig.[BJ we show the maximum value of the Josephson 
current J c for n — 2 and zq = 5, where m is the TKNN 
number on the right superconductor. The two super- 
conductors belong to spin-singlet symmetry. The results 
for m = —4, -2, 0, and 2 in Fig. [UJa) monotonically in- 
creases with decreasing temperature. On the other hand, 
the results for m = 4, 6, and 8 in Fig.^b) show the non- 
monotonic temperature dependence. In particular, the 
results for m = 6 indicate the transition from state to 
7T state with decrease of temperature around T = 0.25T C . 
We also note that the ir state is stable for m = —4 and 8 
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FIG. 6. (color online). The Josephson critical current is plot- 
ted as a function of temperature for zq — 5. We fix TKNN 
number at n = 2 in the left superconductor. The TKNN 
number in the right superconductor is m 

all temperature range. The nonmonotonic behavior of J c 
indicate the instability of the junction between and tt 
states. For m = 2, the Josephson current takes its max- 
imum at <p = 0.55-7T at T = 0. In another cases, CPR is 
almost sinusoidal. The amplitude of J c becomes smaller 
for larger \n — m\ because the integrand of Eq. (1551) oscil- 
lates more frequently as a function of 7 for larger \n — m\. 

In Fig. we show the maximum value of the Joseph- 
son current J c for n = 1 and zq = 5, where to is the 
TKNN number on the right superconductor. The results 
for m = —3, -1, f, and 3 in (a) monotonically increase 
with decreasing temperature. However, the results for 
m = 7 show the nonmonotonic dependence on temper- 
ature because of the 0-tt transition around T = 0.2T C . 
Consistently with the Table. II, the state is stable at 
to =/= 7 and zq 3> 1 all temperature range. 
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